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Abstract 

The  temperature  profile  generated  by  a  thermo-electro-chemical  model  was  used  to  calculate  the  thermal  stress  distribution  in  a  tubular  solid 
oxide  fuel  cell  (SOFC).  The  solid  heat  balances  were  calculated  separately  for  each  layer  of  the  MEA  (membrane  electrode  assembly)  in  order  to 
detect  the  radial  thermal  gradients  more  precisely.  It  appeared  that  the  electrolyte  undergoes  high  tensile  stresses  at  the  ends  of  the  cell  in  limited 
areas  and  that  the  anode  is  submitted  to  moderate  tensile  stresses.  A  simplified  version  of  the  widely  used  Weibull  analysis  was  used  to  calculate 
the  global  probability  of  survival  for  the  assessment  of  the  risks  related  to  both  operating  points  and  load  changes.  The  cell  at  room  temperature 
was  considered  and  revealed  as  critical.  As  a  general  trend,  the  computed  probabilities  of  survival  were  too  low  for  the  typical  requirements  for  a 
commercial  product.  A  sensitivity  analysis  showed  a  strong  influence  of  the  thermal  expansion  mismatch  between  the  layers  of  the  MEA  on  the 
probability  of  survival.  The  lack  of  knowledge  on  mechanical  material  properties  as  well  as  uncertainties  about  the  phenomena  occurring  in  the 
cell  revealed  itself  as  a  limiting  parameter  for  the  simulation  of  thermal  stresses. 

©  2005  Elsevier  B.V.  All  rights  reserved. 
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1.  Introduction 

The  fossil  fuel  resources  mainly  used  for  power  production 
are  limited  and  cause  environmentally  harmful  emissions.  Fur¬ 
thermore,  its  use  is  suspected  to  cause  global  warming.  There¬ 
fore,  its  judicious  use  appears  important.  The  solid  oxide  fuel 
cell  (SOFC)  is  a  direct  conversion  process,  which  allows  the 
production  of  electricity  with  high  efficiency  while  maintaining 
pollutant  emission  at  low  level.  Moreover,  the  SOFC  technol¬ 
ogy  is  now  far  beyond  the  theoretical  state.  Indeed,  prototypes 
have  proven  its  ability  to  achieve  the  conversion  of  fuel  into 
electricity  with  high  efficiency,  especially  in  the  case  of  a  SOFC 
and  gas  turbine  hybrid  system.  The  next  step  is  commercializa¬ 
tion,  which  may  imply  a  drop  of  the  specific  cost  and  extended 
lifetime. 
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The  major  limitation  to  lifetime  comes  from  the  degradation 
of  the  materials,  which  is  due  to  the  aggressive  environment 
related  to  operation  conditions.  A  decrease  of  the  performance 
of  the  cell  is  observed  as  well  as  structural  failure.  The  knowl¬ 
edge  of  the  stress  field  in  the  cell  is  therefore  of  main  interest. 
The  sintering  phase  during  the  manufacturing  process,  com¬ 
bined  with  the  mismatch  between  the  coefficients  of  thermal 
expansion  (CTE)  of  the  materials  of  the  membrane  electrode 
assembly  (MEA)  induces  residual  stresses  in  the  cell.  Dur¬ 
ing  operation,  thermal  gradients  cause  additional  stresses  and 
chemically  induced  stresses  occur  if  interconnectors  sensitive 
to  oxygen  activity  gradient  are  used.  The  consequence  of  these 
phenomena  is  mainly  delamination  and  thermal  cracking  in  the 
critical  layers  of  the  MEA.  Additional  research  is  needed  in  order 
to  quantify  and  predict  the  contribution  of  the  thermal  stresses 
to  the  failure  of  the  cell.  The  operating  condition  of  the  cell  can 
be  adjusted  by  a  suitable  control  strategy.  The  challenge  is  to 
maximize  in  particular  the  efficiency  with  respect  to  lifetime. 

Two  types  of  studies  of  thermal  stresses  in  SOFC  can  be  found 
in  the  literature.  The  first  one  focuses  on  the  residual  stresses  in 
the  cell  at  room  temperature.  Finite-element  analysis  is  com¬ 
pared  to  experimental  measurements  in  order  to  estimate  its 
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Nomenclature 

e  thickness  (m) 

E  Young’s  modulus  (Pa) 

m  Weibull  modulus 

Ps  probability  of  survival 

Pf  probability  of  failure 

r  radius  (m) 

u  displacement  (m) 

Vr  reference  volume  (m3) 

Greek  letters 

a  coefficient  of  thermal  expansion,  CTE  (K-1) 

X  first  Lame  coefficient 

fi  second  Lame  coefficient 

v  Poisson  coefficient 

Oi  principal  stresses,  i=  1,  2,  3  (Pa) 

(j o  characteristic  strength  (Pa) 

r  temperature  difference  (K) 

£2  boundary 

Indices 

i  principal  directions 

j  radial,  axial  direction 

r  radial  direction 

z  axial  direction 


suitability  for  stress  simulation.  Such  studies  were  performed  by 
Yakabe  et  al.  [1],  Atkinson  and  Selguk  [2]  and  Montross  et  al.  [3]. 
The  second  situation  of  interest  is  operating  conditions,  which  is 
known  to  cause  the  failure  of  the  cell.  Lor  example,  Yakabe  et  al. 
[4,5]  and  Selimovic  et  al.  [6]  have  coupled  their  thermo-electro¬ 
chemical  model  with  a  finite-element  tool  in  order  to  simulate 
the  stresses  in  the  cell  during  operation. 

Little  research  has  been  conducted  in  thermal  stresses  in 
SOLC  up  to  now,  compared  to  the  ones  in  thermo-electro¬ 
chemical  modeling.  Lurthermore,  there  are  many  degrees  of 
freedom  in  design,  such  as  tubular  or  planar,  anode  or  cathode 
supported.  Therefore,  no  method  to  handle  the  problem  has  yet 
been  clearly  established  which  is  generic,  reliable  and  afford¬ 
able  in  terms  of  modeling  capabilities.  Moreover,  knowledge  on 
the  mechanical  properties  of  the  materials  of  the  MEA  is  sparse. 
Other  uncertainties  remain  in  critical  areas,  such  as  the  tempera¬ 
ture  at  which  no  internal  stresses  occur  or  behavior  of  the  bonds 
at  the  interfaces  of  the  layers. 

The  purpose  of  the  present  work  is  to  study  the  effects  of 
the  characteristics  of  the  temperature  profile  on  the  stress  field 
in  a  tubular  cell  during  operation.  The  approach  is  to  couple  a 
thermo-electro-chemical  model  for  the  computation  of  the  tem¬ 
perature  profile  with  an  additional  module  for  the  stress  field. 
No  feedback  from  the  stress  calculation  to  the  thermo-electro- 
chemical  model  exists.  Particular  attention  was  paid  to  the  values 
of  the  radial  thermal  gradients  in  the  cell,  which  required  sep¬ 
arate  calculation  of  the  solid  heat  balances  for  the  MEA.  The 
Weibull  analysis  is  a  common  way  to  estimate  the  failure  proba¬ 
bility  of  ceramic  materials.  Several  authors  [3,6]  mentioned  that 


the  local  maxima  of  computed  stresses  might  exceed  the  yield 
strength  of  the  materials  if  realistic  probabilities  of  survival  are 
targeted.  In  this  work,  a  simplified  Weibull  analysis  was  applied, 
which  allowed  studying  the  evolution  of  the  probability  of  sur¬ 
vival  for  several  cases  including  both  steady  state  and  dynamic 
operation.  An  analysis  of  the  sensitivity  of  the  stresses  on  the 
mismatch  between  the  layers’  coefficients  of  thermal  expansion 
was  performed  in  order  to  provide  indication  of  the  mechanical 
requirements  on  the  choice  of  SOLC  materials. 


2.  Modeling 

2.7.  Mechanical  model  and  assumptions 

The  thermo-electro-chemical  SOLC  model,  which  provides 
the  temperature  profile  in  the  cell,  has  been  described  by  Stiller 
et  al.  [7].  The  model  was  implemented  in  gPROMS,  an  equation 
oriented  process  modeling  tool  [8].  It  was  enhanced  in  order  to 
detect  the  radial  thermal  gradients  with  higher  accuracy.  Thus, 
the  MEA  model  is  improved  by  calculating  each  layer  sepa¬ 
rately.  The  steam  reforming,  water  gas  shift  and  electrochemical 
reactions  are  assumed  to  take  place  at  the  interface  between  the 
electrolyte  and  the  anode.  The  model  of  the  tubular  cell  geometry 
is  simplified  to  an  axisymmetric  one,  which  implies  neglect¬ 
ing  circumferential  temperature  variations  and  the  interconnect. 
Lig.  1  depicts  the  geometry  and  provides  the  actual  dimensions. 

The  stress  field  model  is  based  on  the  static  uncoupled  lin¬ 
ear  thermoelastic  theory.  The  modeling  has  been  performed 
with  gPROMS  in  order  have  both  the  thermo-electro-chemical 
model  and  the  stress  model  in  the  same  computational  tool.  This 
approach  is  convenient,  especially  during  dynamic  operation. 
The  drawback  is  the  limitation  of  gPROMS  in  modeling  com¬ 
plex  geometries  and  its  lack  of  advanced  meshing  possibilities. 

The  modeling  assumptions  to  keep  computing  time  reason¬ 
able  are 
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Fig.  1.  Axisymmetric  geometry  of  the  cell  used  for  modeling. 


Table  1 

Material  properties  of  MEA  components 
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Young’s  modulus 
(GPa) 

Poisson’s 

coefficient 

Coefficient  of  thermal 
expansion  (x  10-6  KA1) 

Characteristic 
strength  (MPa) 

Weibull 

modulusa 

Cathode  (LSM) 

35  [20] 

0.25  [20] 

11.7  [17] 

52  [12] 

7  [12] 

Electrolyte  (YSZ) 

1073  K 

183  [12] 

0.313  [12] 

10.8  [17] 

256  [13] 

5  [13] 

298  K 

212  [14] 

0.32  [12] 

332  [13] 

5  [13] 

Anode  (Ni-YSZ) 

57  [14] 

0.28  [12] 

12.2  [17] 

115.2  [14] 

6  [14] 

a  Only  integer  values  could  be  implemented  in  gPROMS. 


The  thermal  expansion  coefficients  are  constant  over  the  con¬ 
sidered  range  of  temperature,  e.g.  values  from  the  literature 
applicable  over  a  range  of  273-1300  K. 

The  values  of  Young’s  modulus  and  Poisson’s  ratio  at  oper¬ 
ating  temperature  are  used  when  available,  but  their  precise 
dependency  on  temperature  is  not  modeled. 

The  value  of  the  zero  stress  temperature1  is  of  the  highest 
importance  for  the  determination  of  the  magnitude  of  the 
stresses.  Major  uncertainties  remain  in  this  area  [1-3].  A  uni¬ 
form  zero-stress  temperature  of  1400  K  is  used  in  the  present 
work. 

The  fixing  of  the  cathode- supported  cell  is  modeled  by  con¬ 
straining  the  axial  displacement  at  the  boundary  Q \  (Fig.  1). 
The  boundary  conditions  at  (Fig-  1)  are  determined  by  the 
pressure  in  the  cell. 

The  seal  (Fig.  1)  at  the  end  is  neglected  due  to  meshing  issues. 
The  boundary  2 3  (Fig.  1)  is  free. 


the  time  being.  They  do  not  necessary  correspond  exactly  to  the 
materials  used  in  tubular  SOFC.  The  input  data  required  for  the 
whole  process  are  listed  in  Table  1. 

2.2.  Verification  against  finite -element  software  (FEM) 

gPROMS  is  not  primarily  designed  to  handle  mechanical 
problems.  Therefore  its  results  were  tested  against  those  pro¬ 
vided  by  a  common  FEM  tool  for  stress  computations,  FEMLAB 
[15].  A  typical  temperature  profile  as  it  appears  during  operation 
was  applied  to  a  tube  of  0.25  m  length,  made  of  layers  of  same 
diameter  and  thickness  as  the  actual  tubular  cell.  Fig.  2  depicts 
the  comparison  of  the  solution  for  the  principal  stresses  of  inter¬ 
est  from  gPROMS  and  FEMLAB.  The  agreement  is  considered 
as  reasonable,  especially  when  considering  that  the  FEMLAB 
model  uses  more  than  40000  nodes,  against  only  738  for  the 
gPROMS  model. 


The  common  Navier’s  equation  is  obtained  by  taking  the  dis¬ 
placement  components  uj  as  the  unknown  of  the  fundamental 
equations  of  linear  elastostatics.  It  is  solved  by  the  model  under 
the  previous  assumptions.  The  resulting  equation  [9,10]  is 


/zV  Uj  +  (A  +  fj)uk,kj 


aE 

- T  7 

1  -  2v  ’} 


The  stress  components  are  computed  from  the  displacement 
fields.  The  three  principal  stresses,2  o\ ,  o 2  and  <73,  are  then 
determined  with  the  principal  invariants3  method  and  used  for  a 
simplified  version  of  the  Weibull  analysis  [11]: 


A  simplification  of  the  choice  of  the  reference  volume,  Vr,  was 
assumed.  It  considers  the  volume  of  the  sample  delimited  by 
the  outer  ring  of  the  so-called  ‘Ting  on  ring  test”,  which  was 
used  in  several  sources  on  material  properties  [12-14].  Only 
few  material  configurations  have  been  mechanically  tested  for 


1  Zero  stress  temperature :  temperature  at  which  no  thermal  stress  act  on  the 
tube.  In  the  SOFC,  this  coincides  with  the  manufacturing  temperature  at  which 
the  layers  are  joined. 

2  Principal  stress :  stresses  in  the  principal  plane  that  is  free  of  shear  stresses. 

3  Principal  invariant,  coefficient  in  the  eigenvalue  problem  required  to  estab¬ 

lish  the  principal  plane,  the  value  of  which  is  not  altered  by  any  rotation. 


2.3.  Investigated  cases 

Both  steady-state  and  dynamic  operation  of  the  SOFC  were 
investigated.  Several  steady-state  cases  were  chosen  from  a 
sequence  of  operating  points  designed  in  studies  on  part-load 
and  load  change  with  a  model  considering  the  SOFC  and  gas 
turbine  hybrid  system  [7] .  These  operating  points  are  character¬ 
ized  by  the  dimensionless  shaft  speed  of  the  gas  turbine  (NDIM). 
High  load  is  denoted  by  values  around  one.  A  particular  feature 
of  the  pursued  operation  strategy  is  that  it  tries  to  keep  the  mean 
temperature  of  the  SOFC  at  a  fairly  constant  value  over  the  range 
of  NDIM  1.00  to  NDIM  0.90.  The  fuel  utilization4  is  close  to 
85%  over  the  whole  range.  Furthermore,  the  stresses  in  the  cell 
at  room  temperature  were  modeled  by  assuming  a  uniform  tem¬ 
perature  of  298  K  in  the  cell. 

The  same  model  is  used  for  dynamic  simulation.  Here,  the 
input  data,  i.e.  operating  conditions,  were  assumed  to  vary  lin¬ 
early  during  a  time  interval  of  60s,  since  the  model  is  limited 
to  the  SOFC  and  does  not  simulate  the  whole  hybrid  system. 
Only  the  SOFC  conditions  corresponding  to  the  load  cases  were 
extracted  from  the  hybrid  system  model.  Two  load  changes  were 
studied,  NDIM  0.65  to  NDIM  1.00  and  NDIM  0.90  to  NDIM 
1.00.  They  represent  changes  of  75  and  15%  of  the  rated  power 
of  the  virtual  hybrid  system,  respectively. 


4 


Fuel  utilization :  fraction  of  the  fuel  which  is  oxidized  in  the  fuel  cell. 
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Gl  FEMLAB 

□  al  gPROMS 

. a2 FEMLAB 

0  a2  gPROMS 

- 03 FEMLAB 

a  a3  gPROMS 

Fig.  2.  Comparison  of  the  thermal  stresses  in  a  tube  computed  by  gPROMS  and 
FEMLAB. 

The  sensitivity  analysis  considers  a  fixed  value  of  the  coeffi¬ 
cient  of  thermal  expansion  of  the  anode,  12.2  x  10  6  K  Awhile 
those  of  the  electrolyte  and  the  cathode  were  allowed  to  vary 
in  agreement  with  data  found  in  the  open  literature.  The  large 
range  of  possible  values  is  due  to  different  doping  methods.  The 
Weibull  modulus  and  strength,  the  Young’s  modulus  and  Pois¬ 
son’s  coefficient  were  assumed  to  be  constant.  The  analysis  was 
performed  at  both  high  and  low  load,  i.e.  NDIM  1.00  and  0.65. 

3.  Results  and  discussion 

3.1.  Stress  fields  during  operation 

A  comparison  of  the  principal  stresses  at  low  and  high  load 
is  depicted  in  Fig.  3.  The  stresses  do  not  vary  significantly  in 


Fuel/airflow  - ► 


Axial  [m] 


- al  at  high  load 

. al  at  low  load 

— b —  a2  at  high  load 

□  a2  at  low  load 

Fig.  3.  Distribution  of  principal  stresses  at  high  and  low  load. 


the  radial  direction.  Two  layers  of  the  ME  A  undergo  tensile 
stresses,  the  electrolyte  and  the  anode.  The  variation  of  the  shear 
stresses  acting  at  the  ends  of  the  cell  induces  high  tensile  radial 
stresses  in  the  electrolyte  at  the  ends  of  the  cell.  They  act  in  very 
small  volume.  The  value  exceeds  200  MPa,  which  is  close  to  the 
characteristic  strength  of  YSZ  at  the  operating  temperature.  This 
behavior  is  opposite  to  that  in  the  middle  of  the  cell,  where  the 
smaller  coefficient  of  thermal  expansion  of  the  electrolyte  leads 
to  compressive  stresses.  In  the  anode,  both  a \  and  02  are  tensile, 
since  the  coefficient  of  thermal  expansion  of  the  Ni-YSZ  cermet 
is  the  highest  among  the  three  layers.  The  stresses  are  highest  at 
the  fuel/air  inlet,  where  the  steam-reforming  reaction  induces  a 
drop  of  the  temperature.5 

3.2.  Steady -state  cases 

The  observation  of  the  specific  values  of  the  stress  field  can¬ 
not  directly  provide  an  indication  of  the  risks  associated  with 
a  particular  operating  point,  since  the  failure  of  SOFC  brittle 
materials  are  related  to  their  defect  structure.  The  probability  of 
survival  is  a  better  criterion.  Fig.  4  depicts  its  value  for  several 
operating  points  of  the  steady-state  sequence.  An  optimum  point 
for  the  global  probability  of  survival  occurs  at  intermediate  load 
and  is  due  to  the  opposite  trends  followed  by  the  probabilities 
of  survival  in  the  layers.  These  are  lower  at  low  load  in  both 


5  In  the  fuel  entering  the  cell,  about  55%  of  the  methane  is  reformed. 
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Fig.  4.  Probabilities  of  survival  as  functions  of  the  operating  point  during  steady  - 
state  operation. 


the  electrolyte  and  the  anode,  due  to  the  higher  temperature 
difference  r  induced  by  the  low  temperature  in  the  cell.  In  the 
cathode,  on  the  other  hand,  the  probability  of  survival  is  lower  at 
high  load,  when  high  radial  thermal  gradients  act,  which  attain 
a  local  maximum  at  NDIM  1.0.  Nevertheless,  it  must  be  noted 


that  the  results  obtained  here  are  restricted  to  a  particular  set  of 
operating  points  and  choice  of  material  properties. 


3.3.  Dynamic  simulations 

The  evolution  of  the  probability  of  survival  in  the  cell  during 
both  large  and  small  load  changes  is  depicted  in  Figs.  5  and  6, 
respectively.  All  load  changes  could  be  divided  into  a  first  phase 
with  quickly  changing  probability  of  survival  and  a  second  phase 
where  the  probability  of  survival  is  slowly  stabilizing. 

From  low  load  to  high  load,  the  increase  of  both  mean  temper¬ 
ature  and  temperature  in  the  reforming  area  at  the  cell  entrance 
induces  an  increase  of  the  probability  of  survival  during  the  first 
phase  of  the  load  change  (Fig.  5).  In  the  second  phase,  the  prob¬ 
ability  of  survival  decreases  due  to  the  eventually  increasing 
thermal  gradients  in  the  reforming  area.  At  the  transition  point 
between  both  phases,  the  conditions  are  still  close  to  the  ones  at 
low  load,  where  the  probability  of  survival  in  the  cathode  does 
not  have  a  major  contribution.  Therefore,  the  trends  are  opposite 
during  the  reverse  process:  the  quick  decrease  of  the  thermal  gra¬ 
dients  is  during  the  first  phase  clearly  effecting  an  increase  of  the 
probability  of  survival.  The  former  is  coupled  to  the  increase  of 
the  temperature  in  the  reforming  area,  which  implies  the  same 
trend  for  the  probability  of  survival  of  the  electrolyte.  In  the 
second  phase,  the  probabilities  of  survival  of  the  anode  and  the 
electrolyte  become  dominant  as  they  follow  the  decrease  of  the 
temperature. 

At  high  load,  the  global  probability  of  survival  is  dominated 
rather  by  the  cathode  than  by  the  anode.  Small  load  changes 
under  these  conditions  are  depicted  in  Fig.  6.  The  influence  of  the 
quick  variations  of  the  thermal  gradients  clearly  appears  on  the 
evolution  of  the  probability  of  survival,  in  both  directions.  The 
influence  of  the  temperature  acts  through  the  same  mechanisms 
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Fig.  5.  Probabilities  of  survival,  temperatures  and  gradients  during  load  changes:  (a)  NIDIM  0.65  to  NDIM  1.00;  (b)  NDIM  1.00  to  0.65. 
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Fig.  6.  Probabilities  of  survival,  temperatures  and  gradients  during  load  changes:  (a)  NIDIM  0.90  to  NDIM  1.00;  (b)  NDIM  1.00  to  0.90. 


as  mentioned  previously.  From  lower  to  higher  load,  the  evolu¬ 
tion  of  the  probability  of  survival  during  the  second  phase  is  less 
dependent  on  the  temperature.  It  is  flatter  during  the  decrease 
of  temperature  from  NDIM  1.0  to  NDIM  0.9.  In  contrast,  the 
evolution  of  the  radial  thermal  gradients  is  almost  symmetric. 
The  previous  results  do  not  show  any  intermediate  drop  of  the 
global  probability  of  survival  during  the  investigated  dynamic 
cases  because  of  the  combined  effects  of  the  temperature  and  the 
radial  thermal  gradients.  This  favorable  combination  may  how¬ 
ever  not  be  generally  true  for  all  imaginable  operation  strategies 
of  the  system. 

3.4.  Sensitivity  on  mismatch  between  coefficients  of  thermal 
expansion 

The  differences  between  the  layers’  coefficient  of  thermal 
expansion  mismatches  are  the  primary  cause  of  thermal  stresses. 
High  uncertainty  exists  about  the  CTE  values.  A  multitude  of 
values  can  be  found  in  the  literature,  because  of  the  wide  possi¬ 
bilities  in  doping.  In  the  present  work,  the  range  of  coefficient 
of  thermal  expansion  was  found  by  literature  review.  In  general, 
the  CTE  may  be  varied  by  modifying  the  material  composi¬ 
tion  or  other  parameters,  such  as  porosity  [16].  Under  these 
modifications,  mechanical  and  electro-chemical  properties  of 
the  materials  would  vary  as  well.  This  dependency  was  how¬ 
ever  not  considered  in  the  present  sensitivity  analysis;  from 
Table  1  only  the  CTE  is  varied  here.  The  CTE  of  the  anode 
was  kept  constant  at  12.2  x  10  6  K  1  [17],  while  those  of  the 
electrolyte  and  the  cathode  were  varied  from  9.7  x  l0-6  K-1 
[17]  to  10.8  x  10-6  K-1  [18]  and  from  11.3  x  lO^K"1  [16]  to 
12.8  x  10-6  K-1  [17],  respectively. 

Fig.  7  depicts  the  results  at  both  high  and  low  load.  The  prob¬ 
ability  of  survival  is  higher  for  high  values  of  the  coefficient  of 


CTE  of  cathode  [lO^K1] 


high  load  NDIM  1.0 
low  load  NDIM  0.65 


Fig.  7.  Effect  of  cathode  and  electrolyte  CTEs  on  the  probability  of  survival  at 
high  and  low  load. 

thermal  expansion  of  the  electrolyte,  since  the  mismatch  relative 
to  both  the  cathode  and  the  anode  is  reduced.  The  optimum  set 
at  low  load  differs  slightly  from  the  one  at  high  load.  Indeed,  the 
cathode  is  the  critical  layer  in  the  latter  case  and  is  prone  to  higher 
tensile  stresses,  if  its  coefficient  of  thermal  expansion  increases. 
In  particular,  the  probability  of  survival  is  critical,  when  the  coef¬ 
ficient  of  thermal  expansion  of  the  cathode  exceeds  that  of  the 
anode.  This  requirement  is  opposite  to  the  one  leading  to  the 
highest  probabilities  of  survival  in  the  anode.  At  low  load,  the 
contribution  of  the  latter  is  larger  and  therefore  induces  a  shift 
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Fig.  8.  Effect  of  cathode  and  electrolyte  CTEs  on  the  probability  of  survival  at 
room  temperature. 

of  the  optimum  point  to  higher  cathode  CTE  compared  with  the 
case  of  high  load.  Furthermore,  highest  probabilities  of  survival 
can  be  achieved  only  for  particular  choices  of  CTE. 

3.5.  Room  temperature  case 

The  difference  between  the  operating  and  the  zero-stress  tem¬ 
perature  has  been  identified  as  a  main  driver  of  thermal  stresses. 
From  this  point  of  view,  room  temperature  appears  as  the  most 
critical  case.  The  shape  of  the  stress  fields  shows  the  same  global 
trend  as  during  operation;  that  is  high  stresses  at  the  ends  of  the 
cell  in  the  electrolyte  and  tensile  stresses  in  the  whole  anode. 

Fig.  8  displays  the  result  of  a  sensitivity  analysis  of  the  same 
kind  as  in  Section  3.4.  The  optimum  probability  of  survival  over 
the  considered  CTE  range  is  very  low  and  located  in  a  narrow 
regime. 

3.6.  Limitations 

The  comparison  of  the  results  presented  in  Sections  3.2-3 .4 
shows  to  some  extent  contradiction  with  the  common  observa¬ 
tions,  which  indicate  the  failure  of  the  cell  during  operation. 
However,  the  computation  of  stresses  performed  by  Montross 
et  al.  [3]  in  the  case  of  a  planar  cell  at  room  temperature  also 
indicates  that  probabilities  of  failure6  as  low  as  10-5  cannot  be 
achieved  with  the  common  SOFC  materials.  Nevertheless,  this 
requirement  is  typical  for  applications  targeted  by  the  SOFC 
technology.  Commercialization  is  obviously  not  possible  under 
these  conditions.  Atkinson  and  Sefimk  [12]  pointed  out  the  low 
values  of  the  Weibull  moduli  as  a  limiting  parameter.  Ideally, 
these  should  be  higher  than  10. 

Mechanical  data  for  SOFC  materials  is  still  sparse.  In  partic¬ 
ular,  the  dependency  of  the  Young’s  modulus  and  the  strength  on 
temperature  has  not  yet  been  precisely  quantified.  Its  exploration 


6  Probability  of  failure:  Pf  =  1  —  Ps. 


Fig.  9.  Weibull  plot  of  75%  Ni(0)-YSZ  cermet  (anode  material)  at  room  tem¬ 
perature:  (a)  unreduced  [12];  (b)  unreduced,  porosity  0.06  [14];  (c)  reduced, 
porosity  0.27  (used  for  the  calculation  in  the  present  work)  [14];  (d)  reduced, 
porosity  0.36  [14];  (e)  reduced,  porosity  0.39  [14]. 

is  complicated  by  the  wide  possibilities  in  material  composition, 
which  mainly  has  to  fulfill  electro-chemical  requirements.  The 
manufacturing  process  also  plays  a  major  role;  suitable  sintering 
temperatures  and  processes  are  under  investigation.  The  exper¬ 
imental  data  therefore  depend  on  the  supplier  of  the  sample. 
The  reduction  of  the  anode  material  is  another  issue.  Only  few 
researchers  [14]  have  consistently  published  its  influence  on  the 
strength.  This  fact  is  visible  in  Fig.  9.  It  shows  the  probability  of 
survival  of  75%  Ni(0)-YSZ  over  the  stress  at  room  temperature. 
The  large  discrepancies  are  partially  due  to  different  porosity  and 
state,  i.e.  reduced  and  unreduced.  However,  in  the  most  compa¬ 
rable  case,  the  strength  of  the  tested  material  differs  significantly. 
For  the  preceding  calculations,  material  (c)  is  used.  The  results 
presented  in  Sections  3. 2-3. 5  obviously  change  if  material  (a) 
or  (b)  was  used.  The  dependency  on  porosity,  which  increases 
after  reduction,  amplifies  the  problem. 

Other  fundamental  uncertainties  remain.  In  particular, 
extended  knowledge  of  the  behavior  of  the  bonds  between  the 
layers,  as  well  as  on  the  phenomena  determining  the  value  of  the 
zero  stress  temperature  are  required.  In  the  case  of  a  tubular  cell 
manufactured  by  electrochemical  vapor  deposition,  the  process 
requires  the  insertion  of  additional  layers  in  order  to  prevent  the 
growth  of  undesired  phases  [17].  A  model  assuming  three  lay¬ 
ers  joined  to  each  other  at  a  uniform  temperature  may  not  be 
refined  enough  to  simulate  the  stresses  in  such  tubular  SOFC 
with  satisfactory  accuracy.  The  work  of  Fi  et  al.  [19]  observed 
complex  behavior  of  cell  samples  during  their  experiment  for 
the  analysis  of  crack  density.  In  particular,  reproducibility  could 
not  be  achieved.  The  effect  of  the  sintering  temperature  on  the 
crack  density  was  clearly  observed.  Unexpected  behavior  of  the 
layers,  such  as  compliance  of  the  electrolyte  is  reported. 

The  experimentally  observed  mechanical  failure  of  SOFC 
depends  on  the  degradation  of  the  materials  due  to  the  aggres¬ 
sive  environment  during  operation.  Studies  on  the  evolution  of 
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the  mechanical  properties  of  the  materials  during  the  life  of  the 
SOFC  are  necessary  to  predict  the  SOFC  life.  Such  phenomena 
have  not  been  considered  in  the  present  study.  The  probabil¬ 
ity  of  survival,  which  has  been  used  as  test  criterion,  considers 
only  the  initial  structure.  Indeed,  the  current  knowledge  on  the 
mechanical  phenomena  in  SOFC  does  not  allow  modeling  for 
the  purpose  of  life  assessment. 

4.  Conclusions 

A  model  for  calculating  stresses  was  added  to  a  thermo- 
electro-chemical  model  in  order  to  simulate  the  thermal  stresses 
in  a  tubular  cell.  For  all  studied  cases,  the  stress  fields  display 
high  tensile  stresses  in  limited  areas  at  the  ends  of  the  cell  in  the 
electrolyte  and  significant  values  in  the  whole  anode.  The  latter 
is  due  to  its  high  coefficient  of  thermal  expansion.  The  mag¬ 
nitude  of  the  stresses  is  higher  at  the  fuel/air  inlet,  where  the 
internal  steam-reforming  reaction  induces  a  temperature  drop. 
The  comparison  with  FEMLAB,  a  suitable  FEM  tool,  confirmed 
the  ability  of  the  modeling  tool  gPROMS  to  handle  mechanical 
problems  in  particular  cases,  despite  its  limitation  in  meshing 
abilities,  i.e.  uniform,  rectangular  meshing  and  simple  geome¬ 
tries. 

The  assessment  of  the  risks  related  to  particular  operating 
points  and  load  changes  was  performed  through  the  probabil¬ 
ity  of  survival.  The  evolution  of  the  latter  is  dependent  on  the 
temperature  difference  in  both  the  anode  and  the  electrolyte, 
whereas  in  the  cathode,  radial  thermal  gradients  are  most  impor¬ 
tant.  Therefore  a  maximum  point  for  the  value  of  the  probability 
of  survival  of  the  cell  is  observed  at  intermediate  load.  Transient 
operation  did  not  lead  to  any  drop  in  the  value  of  probability  of 
survival.  However,  this  favorable  behavior  may  change  when  a 
different  operation  strategy  is  applied. 

The  greatest  difference  between  the  actual  temperature  and 
the  sintering  temperature  occurs  when  the  cell  is  at  room  tem¬ 
perature.  The  probability  of  survival  calculated  here  reaches 
critical  values,  which  might  be  acceptable  for  prototypes,  but 
does  not  allow  any  commercialization.  Significant  improvement 
of  the  mechanical  properties  of  the  materials  is  required.  This 
result  is  somewhat  contradictory  to  the  common  observations, 
which  rather  indicate  failure  during  operation.  As  a  matter  of 
fact,  it  shows  the  limitation  of  the  thermal  analysis  alone  and 
points  out  the  importance  of  considering  the  degradation  of 
the  material  during  operation.  Unfortunately,  it  appears  that 
the  published  knowledge  on  material  properties  and  complex 
phenomena  occurring  in  the  cell  is  not  sufficient  to  allow  any 
modeling  for  the  purpose  of  lifetime  assessment.  Therefore,  fur¬ 
ther  studies  should  focus  on  the  initiation  and  growth  of  cracks 
in  SOFC  and  on  improving  the  knowledge  of  material  proper¬ 


ties.  An  in-depth  study  of  the  influence  of  the  choice  of  control 
strategy  would  also  be  of  interest  in  order  to  find  general  criteria 
for  stress-friendly  operation. 
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